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Abstract. In this paper, we propose a probabilistic approach to the 
study of the characteristic polynomial of a random unitary matrix. We 
recover the Mellin Fourier transform of such a random polynomial, first 
obtained by Keating and Snaith in [7], using a simple recursion formula, 
and from there we are able to obtain the joint law of its radial and an- 
gular parts in the complex plane. In particular, we show that the real 
and imaginary parts of the logarithm of the characteristic polynomial 
of a random unitary matrix can be represented in law as the sum of in- 
dependent random variables. From such representations, the celebrated 
limit theorem obtained by Keating and Snaith in 7 is now obtained 
from the classical central limit theorems of Probability Theory, as well 
as some new estimates for the rate of convergence and law of the iterated 
logarithm type results. 



1. Introduction 

In [7], Keating and Snaith argued that the Riemann zeta function on the 
critical line could be modelled by the characteristic polynomial of a random 
unitary matrix considered on the unit circle. In their development of the 
model they showed, via calculating the Mellin-Fourier transform, that the 
logarithm of the characteristic polynomial weakly converges to a normal dis- 
tribution, analogous to Selberg's result on the normal distribution of values 
of the logarithm of the Riemann zeta function |17j . 

The purpose of this paper is to prove an equality in law between the char- 
acteristic polynomial and products of independent random variables. Using 
this we rederive the limit theorem and Mellin-Fourier transform of Keating 
and Snaith and prove some new results about the speed of convergence. 

Let Vjv denote a generic N x N random matrix drawn from the unitary 
group U (N) with the Haar measure ^t/(7v) . The characteristic polynomial 
of Vn is 

Z{Vn,0) := det(Ijv - e-'Vjv) 

N 



n 
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where e^^^ , • • • , e'^^ are the eigenvalues of Va^ • Note that by the rotation 

invariance of Haar measure, if 6 is real then ZiV^, 9) Z{yN,Q). Therefore 
here and in the following we may simply write for Z{Vn,9). Keating 
and Snaith [7] evaluated the Mellin- Fourier transform of Z^. Integrating 
against the Weyl density for Haar measure on U{N), and using certain 
Selberg integrals, they obtained, for all t and s with DKe{t is) > —1, 

In ItTj and in this article, arg Z^ is defined as the imaginary part of 

N 

log Zn := J^log(l-e*^") 

n=l 

with Jmlog(l — e'^") G (— 7r/2, 7r/2]. An equivalent definition for log^Ar is 
the value at point x = 1 of the unique continuous function log det{lN — xVjy) 
(on [0, 1]) which is at x = 0. 

By calculating the asymptotics of the cumulants of (jl.ip , they were able 
to show that for any fixed s,t, 

E (^\ZN\'/V(^^^^'e''''<^^^/V^^^^) ^ exp {^t^ - ^s^) 

as — > oo, and from this deduce the central limit theorem 

log Zn law 



ilogA 



AAi + i7V2, 



where A/i and J\f2 are two independent standard Gaussian random variables. 

We will see in the two following sections how (1.1) may be simply in- 
terpreted as an identity in law involving a certain product of independent 
random variables. In particular, we shall show that 9^e log Zjv and 3m log Z^ 
can be written in law as sums of independent random variables. Sums of 
independent random variables are very well known and well studied objects 
in Probability Theory, and we can thus have a better understanding of the 
distribution of the characteristic polynomial with such a representation. We 
also apply the classical limit theorems on such sums to obtain asymptotic 
properties of Z^ when N — > oo. In particular, we recover the convergence 

in law of log ZAr/y'^^dogiV to a standard complex Gaussian law as a con- 
sequence of the classical central limit theorem. We also obtain some new 
results about the rate of convergence and prove an iterated logarithm law 
for the characteristic polynomial. 

More precisely, the paper is organized as follows: in Section [21 we use a 
recursive construction for the Haar measure on U (N), to obtain our first 
equality in law for the distribution of the characteristic polynomial as a prod- 
uct of independent random variables, from which we obtain a new proof of 
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(jl.ip which does not use Selberg's integrals or the Weyl density. Then in Sec- 
tion[3]we use (jl.ip to deduce the joint law of (9^e log Zat, Jmlog Zn), writing 
each component as a sums of independent random variables. Using these 
two representations for Z]\[, in Section |4] two new proofs of Keating-Snaith 

limit theorem (the convergence in law of log Z7v/Y^^^"^og~/V to a standard 

complex Gaussian law) are provided. We also give estimates on the rate of 
convergence in the central limit theorem. In Section [5] we see how strong 
limit theorems such as the iterated logarithm can be deduced from our rep- 
resentations. 

In a companion paper to [7j, Keating and Snaith [8] studied the char- 
acteristic polynomial for classical compact groups other than the unitary 
group, and in Section [6l we also give similar results for SO{2N) which plays 
a similar role to U{N) for other families of L-functions. 

Since the publication of [7j, there have been many developments in un- 
derstanding the distribution of the characteristic polynomial. For example, 
other limit theorems and large deviation results were derived for the char- 
acteristic polynomial by Hughes, Keating and O'Connell in [5|. The distri- 
bution of the characteristic polynomial away from the point = in the 
other groups has been studied by Odgers [12]. For more details about the 
connections between random matrix theory and analytic number theory, see 
|10j and the references therein, or the excellent survey article by Royer [16j . 

2. Decomposition of the Haar measure 

In this section, we give an alternative proof of formula (jl.ip . which does 
not necessitate the explicit knowledge of the Weyl density formula and the 
values of some Selberg integrals. This new demonstration relies on a recur- 
sive presentation of the Haar measure ^iu{N)- 

2.1. Decomposition of the Haar measure. Let Vn be distributed with 
Haar measure Hu{N) on U{N). If M G U{N + 1) is independent of Vn a 
natural question to ask is under which condition on the distribution of M, 
is the matrix 

" {I I) P-'> 

distributed with the Haar measure in dimension + 1 ? 

The solution to such a question allows one to recursively build up a Haar 
distributed element of U{N). 

This question was partially answered by Mezzadri relying on a general 
method due to Diaconis and Shahshahani [3]. Due to a factorization of the 
Ginibre ensemble, Mezzadri showed that when M is a suitable Householder 
reflection then (j2.ip is distributed with the Haar measure in dimension A^-|-l. 
More precisely, he showed that if t; is a unit vector chosen uniformly on the 
{N + l)-dimensional unit complex-sphere 
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and if 9 is the argument of the first coordinate of v, and u is the unit vector 
along the bisector of ei and where ei = (1, 0, . . . , 0) is the unit vector for 
the first coordinate, then one could take M to be an element of U{N + 1) 
that can be written — e~'^(/jv+i — 2uu^). 

By the application M i— *■ Mi, it is clear that a necessary condition for 
our question must be that v must be distributed according to the uniform 
measure on S^^^"^ . The following proposition states that this condition is 
also sufficient. It is a slight generalization of Mezzadri's result, and its proof 
does not require a decomposition of the Ginibre ensemble. 

Proposition 2.1. Let Me U{N + 1) be such that its first column Mi is 
uniformly distributed on 5^^'^'^ . IfVjsf G U{N) is chosen independently of 
M according to the Haar measure fJ-jjiN)) then the matrix 

is distributed with the Haar measure tJ-u{N+i) ■ 

Proof. Due to the uniqueness property of the Haar measure, we only need 
to show that for a fixed U G U{N + 1) 

In the following, a matrix A will often be written where Ai is its 

first column. As ?7 € U{N + 1), (C/M)i = UMi is distributed uniformly on 
the complex unit sphere y^^'^, so we can write UM = (Pi ||P), with Pi uni- 
formly distributed on and P having a distribution on the orthogonal 
hyperplane of Pi . We then need to show that 

where all matrices are still independent. As Mi and Pi are identically dis- 
tributed, by conditioning on Mi = Pi = -y (here v is any fixed element of 
•^C^^) is sufficient to show that 

(.IIP') {ID (.||M') {II), 

where M' (resp P') is distributed like M (resp P) conditionally to Mi = 
V (resp Pi = v). Let A be any element of U{N + 1) such that A{v) = 
(1,0,..., 0). Since A is invertible, we just need to show that 

A(vnp'){i vl)=-^Hi"')(; 

that is to say 

P"Vn '= M"Vn, 
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where P" and M" are distributed on U{N) independently of Vn- By in- 
dependence and conditioning on P" (resp M"), we get P"Vn ^= Vn (resp 



The result of this proposition is very natural. It states, roughly speaking, 
that in order to choose uniformly an element of U {N + 1) (that is to say an 
orthogonal unitary basis) one just needs to choose the first element uniformly 
on the sphere and then an element of U{N) in the hyperplane orthogonal 
to the first element, uniformly. 

2.2. Recovering the Mellin Fourier transform. The decomposition of 
the Haar measure presented in the previous paragraph gives another proof 
for equation (jl.ip . In reality, the following Proposition [22] gives much more, 
as we get a representation of Z{Vn) as a product of N simple independent 
random variables. 

Proposition 2.2. Let Vn G U{N) be distributed with the Haar measure 
I^U{N) ■ Then for all 9 £M. 



with 6*1, ... , On, /3i,0) • • • ! Pi,n-i independent random variables, the Ok 's uni- 
formly distributed on [0, 2tt] and the Pij 's (0 < j < N — 1) being beta 
distributed with parameters 1 and j (by convention, (5ifi is the Dirac distri- 
bution on 1). 

Proof. As previously mentioned, it suffices to consider the case = 0. 

Note now that in Proposition l2.lt we can choose any matrix M G U{N) 
with Ml uniformly distributed on the complex sphere ■ Let us choose 
the simplest suitable transformation M : the reflection with respect to the 
median hyperplane of ei and Mi, where Mi is chosen uniformly on ■ 
Let the vector v be Mi — ei. Therefore there exists (A2, . . . ,\n) G C^~^ 
such that 



M"Vn = Vn) by definition of the Haar measure ij.u{N)- This gives the 
desired result. □ 



N 



det(I^ - c^'Vn) n (1 + e''^ 



k=l 



M = {ei + v\\e2 + \2v\\ 
So, with Proposition 12.11 one can write 



CN + Xnv) . 




M 
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If we call . . . ||ti7v-i) := ^iv_i then using the multi-linearity of the 
determinant we get 



det 



1 

FLi 



M 



det I —v\ 
det 
det 





Ul 



Ul 



62 - X2V\\ . . . 





UN-l 



62 . . 





UN-l 



-Vl 



— T 



Finally, 



-Vl det - In-1 



det{lN - Vn) ^= -Vl det(/Ar_i - Va^-i 



with —Vl = 1 — Mil and Vn-i independent. Therefore, to prove Proposi- 
tion 12.21 we only need to show that Mu ^= e'^^ y^^VV-i"- This is straight- 
forward because, since Mi is a random vector chosen uniformly on , we 
know that 



law 



11 



xi + iyi 



JSn fp 



with the Xi^s and y^'s all independent standard normal variables, On and 
Pi^N-i as stated in Proposition 12. 2i □ 

To end the proof of (jl.ip . we now only need the following lemma. 

Lemma 2.3. Let X := 1 -\- e'^y^, where 6 has uniform distribution on 
[0, 27r] and, independently (3 has a beta law with parameters 1 and N — \. 
Then, for all t and s with 9^e(t is) > —1 

r (iV) T{N + t) 



E {\X 



t isnvgX 



) 



r (AT + *±^) r (TV + 

Proof. First, note that 

E {\X\'e''^'^^) = E (^x'^t+s)/2Y(*-'y^'^ 

E( {l + e''^y (l + e 



-ie 



with a = {t + s)/2 and b = (t — s)/2. Recall that if |x| < 1 and u G M then 



{i+xr = Y, 



k=0 



u{u — 1) . . . {u — k + 1) ^ 

k'. 



E 

A;=0 



kl ^ ' 
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where {y)k = y{y+l) . . . (y+k — l) is the Pochhammer symbol. As \e^^\f/3\ < 
1 a.s., we get 



E 



kl 



-b)e r,L _ 



After an expansion of this double sum (it is absolutely convergent because 
*>le(t ± s) > —1), all terms with k ^ t will give an expectation equal to 0, 
because we integrate with respect to the uniform probability measure along 
the unit circle. So we get 



E 



{-a)k{-h)k ok 



.A;=0 



(k\ 



As /3 is a beta variable with parameters 1 and — 1, we have 



E 



r(l + A;)r(A^) k\ 



hence 



r(l)r(Ar + A;) (Ar)fc' 
°° (-a)fc(-6)fc 



fc=o 



k\{N), 



Note that this series is equal to the value at z = 1 of the hyper geometric 
function H(—a,—b,N;z). This value is well known (see, for example, [1]) 
and yields: 

n Yitjs.rsxi _ r{N)T{N + a + b) 

' J r(A + a)r(A + 6)- 

This gives the desired result. □ 

Comments about Selberg integrals. To prove (jl.ip Keating and Snaith [7J , 
relying on Weyl's integration formula, used the result by Selberg 



J(a,6,a,/3,7, A) := 



TV 



(2tt 



(a + 5)(a+/3)Ar-7Af(TV-l)-Af 



N-1 

n 

j=0 



r(l + 7 + j7)r(« + /3( A + j - 1)7 - 1) 
r(l+7)r(a-i7)r(/3-j7) 



(2.2) 

in the specific case a = b = 1 and 7=1. Thus our probabilistic proof 
for (II. ip also gives an alternative proof of (j2.2p for these specific values 
of the parameters. Moreover, as we will see in Section [6l a similar result 
as Proposition 12.21 still holds for the orthogonal group. As a consequence 
Selberg's integral formula admits a probabilistic proof for a = 6 and 7=1/2 
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and 1 (however this method rehes on Weyl's integration formula, which is 
essentially analytical) . 

2.3. Decomposition of the characteristic polynomial off the unit 
circle. Proposition 12.21 can be extended to the law of the characteristic 
polynomial of a random unitary matrix off the unit circle where we replace 
e*^ by a fixed x. Once more, due to the rotation invariance of the unitary 
group, we may take x to be real. 

Proposition 2.4. Let € M, Vn~i distributed with the Haar measure 
/i[/(Ar_i), Ml uniformly chosen on -5^^ , independently of Un-i. We write 
Mil for the first coordinate of Mi, and Mi for the vector with coordinates 

Mi2,...,MiN. 

Then, ifV^ is distributed with the Haar measure fJ^u{N)> 



det(/ 



N 



xV^) (1 



x(l — x)-~-T —- 



xMii) 

■T 



-1 



det(/Ar„i - xVn~i) 
- xIn-i)^^Mi det (/at- 1 



xVn-i). (2.3) 



1-Mn 

Proof. The idea is the same as for Proposition 12.21 where we use Propo- 
sition 12.11 with a specific choice of M, the reflection with respect to the 
hyperplane median to ei := (1,0, . . . , 0) and Mi, a vector of chosen 
uniformly. If k := Mi — ei we can write more precisely 



M 



I -Mil I -Mil 



Thus, using multi-linearity of the determinant, due to Proposition 12.11 we 
get after some straightforward calculation 



det (I 



N 



xVn) ^= det 



b det 




with b = and a = xMii){i Mn) ^ ^ want to express these 

terms with respect to det(/Ar_i — Vn-i), writing B := Vj^_i — x/jv-i leads 
to 

detd. - XV.) b det ( f ) det ( ) 

= bdet( a-MiB-^Mi ••• ^ 
V BVn-i I 

= 6 (a - Ml B-^Mi) det(/Ar_i - xVn-i). 
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This is the expected result. □ 

One may try to get a new proof of Weyl's integration formula thanks to 
this recursive construction of a characteristic polynomial. Let un be the 
probability measure on [0, 27r)^ with density 

ly^idai, . . . ,daN) = cat je""^ - e*"'=|^dai . ..duN- 

It would be sufficient to show that if {6i, . . . ,9n) and (^i, . . . , ^at-i) are 
independent and respectively distributed according to i^j\f and i'n-i, then 
for all X G R, with the notations of the proposition, 

N N-l 



However, this identity in law does not seem to have an easy direct explana- 
tion. 

3. Decomposition into independent random variables 

3.1. Some formulae about the beta-gamma algebra. We recall here 
some well known facts about the beta-gamma algebra which we shall often 
use in the sequel. A gamma random variable 7^ with coefficient o > has 
density given by: 

P{7a G dt} = ^e-'dt. 
T{a) 

Its Mellin transform is (s > 0) 

^ r{a + s) 

"'^J r(a) 

A beta random variable fia,b with strictly positive coefficients a and b has 
density on [0, 1] given by 

^{Pa,,edt} = ^j^t'^-^{l-tr'dt. 

r (a) r (6) 

Its Mellin transform is {s > 0) 

^ ^^^''^ - r(g) r(g + 6 + .)- ^^-^^ 

We will also make use of the following two properties (see |3] for justifi- 
cations): the algebra property where all variables are independent 

^ law 
Pa,b7a+b = la, 
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and the duplication formula for the gamma variables, with all variables 
independent 



law „ / / 

7i = 2j7i7i±i- 



3.2. The joint law of (|.^Ar|, CJmlog Zj^). In this Section, we use the Mellin 
Fourier transform (jl.ip obtained in Section [2] to deduce the joint law of 
(|ZAr|, 3mlog Ztv). For simplicity, let us write 

Atv = I Zjy I , and Ijy = 3m log Zjy 

so with this notation formula (jl.ip states 

Lemma 3.1. Let Wj have density 
where 

220-1) 



and let 



^(2i-2)! 



X := f3jj^i2cosWje^i. 



where all the random variables in sight are independent. Then 

LI I J r{j + {t + s)/2)T{j + {t-s)/2) 

Proof. By the definition of X, we have that 

E [iXfe'''-''^^] = 

f7r/2 



(3.3) 



J-tt/2 

By p.ip we have 



r(i)r(2i - 1 + 1) 

Note that 



= (l + g2ix-)i-l+{*+s)/2 _^ g-2ix^i-l+(t-s)/2 

Both terms on the RHS can be expanded as a series in e^'^ or e^^'^ for all 
X other than x = 0. Integrating over x between — 7r/2 and it/2, only the 
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diagonal terms survive, and so 

r.7r/2 
'-7r/2 

{-{j - 1 + (t + s)/2)), (-(j - 1 + (t - s)/2)). 



J-7r/2 



k\k\ 
k=0 

= H{-{j -l + {t + s)/2), -(j - 1 + (t - s)/2), 1; 1) 

where H is a hypergeometric function. The values of hypergeometric func- 
tions at z = 1 are well known (see, for example py), and are given by 

r(2j - 1 + 1) 

r(j + {t + s)/2)r{j + {t- s)/2) 

and this completes the proof. □ 

The next Theorem now follows easily from the previous lemma: 

Theorem 3.2. Let An = \Zn\ o^nd 1^ = '3mlogZj\[. Let (/?jj-i)i<j<^ 
be independent beta variables of parameters j and j — 1 respectively (with 
the convention that /3i_o = Ij- Define Wi, . . . ,Wj\[ as independent random 
variables which are independent of the {f3jj-i)^^j^j^, with Wj having the 
density: 

^2U-i) (dv) = ''^V(2f-"2)?'^' (3.4) 
Then, the joint distribution of (/at, A^r) is given by: 

(N N \ 

i=i i=i / 



We now recover a formula obtained in [TT] in the study of the relations 
between the Barnes function and generalized gamma variables. To this end, 
we need the following elementary lemma whose proof is left to the reader: 

Lemma 3.3. Let Vt be distributed as 

F{Vtedv) = ^ \v\<-, 

with Ki = . Then 

(r(i+|))' 

cos {Vt) ^1 0t+i 1 . 

IfWj '= V2{j-i) then 
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Proposition 3.4 ([TT]). Let (7j)i<j<7v ^''^^ v-? j sequences of in- 

dependent gamma(j) variables. Then we have 

N N 
law 



117.'= A^nV^n^ (3.7) 

Proof. Considering only the second component in (j3.5p . and multiplying 

N 

i= 



both sides by Hy^i 72j-i) ^^'^ thanks to the beta-gamma algebra, we obtain: 



A., 72.-1 j = [U ij 2^ [U cos iW,)j . (3.8) 

Now we apply the lemma to the right hand side of the above equality in law 
to obtain that 

N \ / N \ /^ \ ^ 

n 7, 2" n - f"'.) = n 7. n ^'■'> 

On the other hand, from the duplication formula for the gamma function, 
we have for any a > 



law 



7a = 2. haj',+i 



^ 2 



thus, on the left hand side of (jS.Sp we get 



N N 
■ law . jv 



A.2- n ./w^ '= A.2- n ^M:7^. (3.10 



2 ' 2 



Now comparing (j3.9p and (j3.10p we obtain 



Af Af 
law 



□ 



Infinitely divisible laws form a very remarkable and well studied family 
of laws in Probability Theory. It is easily see from Proposition 13.41 that the 
law of log \ Zn\ is infinitely divisible. 

Proposition 3.5. The law o/log|Zjv| is infinitely divisible. 

Proof. It follows from Proposition 13.41 and the fact that logarithm of a 
gamma variable is infinitely divisible (see, for example, [2]). □ 

Remark. The law of Jm log Z]\f is not infinitely divisible since it is a bounded 
random variable. 
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4. Central limit theorems 

In this section, we give two alternative proofs of the following central limit 
theorem by Keating and Snaith [7]. The first one from the decomposition in 
Section [21 the second from the last decomposition in Section [3l The original 
proof by Keating and Snaith relies on an expansion of formula (jl.ip with 
cumulants. 

Theorem 4.1. Let '■= det(/Ar — Vat), where Vn is distributed with the 
Haar measure on the unitary group U (N) . Then, 

^?i£^^Af, + iAf2, (4.1) 
i log TV 

as N ^ oo, with A/i and J\f2 independent standard normal variables. 

4.1. Proof from the decomposition in Section [2], From Proposition l2.2t 
we know that 

N 

dei{lN - Vn) '= n (l + e''^- v^^^) , 

k=l 

with 6*1, ... , On, Pi,o, ■ ■ ■ , (3i,n-i independent random variables, the 9k s uni- 
formly distributed on [0,27r] and the /3i,j's (0 < j < — 1) being beta 
distributed with parameters 1 and j. 
In the following we note 

N 

Xn := J2 log (l + e''' VK^) > 

i=l 

with log(l + e) defined here as X^j>i(— l)'''''^e-'/i (this is convergent a.s. 
because \e^^'= (3i^k-i\ < 1 a.s.). 

What we need in order to prove Theorem 14.11 is the following : 

(1) first to show that is equal in law to logZjv, as it is defined in 
the introduction. It is not so obvious, because the imaginary parts 
could have a 2kTr difference. 

(2) then show that Xn converges to what is expected in (4.1). 

Proof for (1). Equation ()2.3p . stated for a fixed x, is also obviously also 
true for a trajectory, for example for x G [0, 1], 

(det(/Ar -xVAr),0 <x<l) 

'= ((l-/(x,y7v-i,Mi))det(/^_i-xyiv-i),0<x< 1), 

with the suitable / from (|2.3|) . Let the logarithm be defined as in the 
Introduction (ie by continuity from x = 0). The previous equation then 
implies, as / is continuous in x, 

logdet(/7V - xVn) '= log(l - /(x,T/^_i,Mi)) + logdet(/,v-i - xVn-i). 
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One can easily check that \ f(x, Vn~i, Mi)\ < 1 for all x G [0, 1] a.s., so 

iog(i - /(x,y^_„M,)) = 

— 7 

for all X G [0, 1] almost surely. In particular, for a; = 1, we get 



logdet(/^^ - Vn) '= \^J.-^y^'^j + logdet(lA,_i - Vn-i), 

which gives the expected result by an immediate induction. We have there- 
fore shown that log Zjsf = Xjv . 

Proof for (2). The idea is basically that /3i,A;-i tends in law to a Dirac dis- 
tribution on as /c tends to oo. So log(H-e'^'= ^f^i,k-i) is well approximated 
by 6'^*= ^J|3l^k-l^ and as this has a distribution invariant by rotation, the cen- 
tral limit theorem will be easily proven from classical results in dimension 
1. 

More precisely, Xjv can be decomposed as 

N N N , .j+i 

fc=l A;=l j>Zk=\ •' 

Xi(iV) X2(Ar) XsiN) 

where all the terms are absolutely convergent. We study these three terms 
separately. 

Clearly Xi{N) has a distribution which is invariant by rotation, so to 
prove that ^'^^^^ A/i + iA/2, we only need to prove the following result 

log AT 

for the real part : 

J2k=i c os 9k-\/Pi,k-i law 

where Af is a standard normal variable. As E(cos^(^fe)/3i,fe-i) = this is a 
direct consequence of the central limit theorem (our random variables check 
the Lyapunov condition). 

To deal with X2{N), as ^jt>ol/^^ < co, there exists a constant c > 
such as E(|X2(Af)|2) < c for all N e N. Thus (X2(A^),iV > 1) is a L^- 
bounded martingale, so it converges almost surely. Hence 



X2{N)/^llogN ^0 a.s. 

Finally, for X^iN), let Y := Ej^ls Er=i One can easily 

check that K{Y) < cxd, so F < 00 a.s., so as AT — 00 

l^3(iV)|/ 5 log < Y/^^ log AT ^ a.s. 
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Gathering all these convergences, we get the expected result : 

law 



i log TV 

with the notations of Theorem 4.1. 



4.2. Proof from the decomposition in section [3l We shall give here 
a very simple proof of the central limit theorem for based on the 

decomposition into sums of independent random variables and a classical 
version of the multidimensional central limit theorem. 
From Theorem 13. 2^ we have: 

pmiogz^,iog|z^l)^^" I yw^yy, i . (4.2) 




where 

Tj = log(/5jj_i2 cos Wj). (4.3) 

Now, from the discussion preceding Theorem 13.21 we have for s > — 1 and 
t > -1: 

E\e''^i] = (4 4) 

^ ^ r(, + i)r(,-f)' ^'-'^ 

From these Fourier transforms, one can easily deduce the moments or the 
cumulants of all orders for Wj and Tj (see \15\ ?, ?] for definition of cumu- 
lants and their relations with moments) by taking successive derivatives at 
0. For our purpose, we will only need the first three moments or cumulants. 
Since the calculation of the derivatives have already been done in [7], we 
will only recap them here. Let us call Qj^k the A;-th cumulant of Tj and Rj^k 
the k^^ cumulant of Wj. Then we have 



and 



where 



R 



if k is odd 
IS even 

,(,),^, _ d'=+Mogr(z) 




are the polygamma functions. Now, since the cumulants of a sum of in- 
dependent random variables are the sum of the cumulants, we can easily 
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obtain that the cumulants of ^jLi and ^^Li ^''^^ respectively 
and 

I if fc is odd 

^fc/2+l 

Moreover, we have the following expansion of the polygamma function (see, 
for example Corollary 1.4.5 of [Ij): 



m-^Ef=i^('-'Mj) if A: is even 



oo „ 

^(.)~logz---^^ (4.6) 

and 



2z ^ 2nz2n 

n=l 



^p(>'){z) = (-1) 



k-l 



{k-iy. kl ^ (2n + A;-l)! 

n=0 



(4.7) 



for \z\ — > oo and | argz| < vr, and where the i?2n are the Bernoulli numbers. 
We deduce from (j4.6p that the variances of 9\t log Z]\f and 3m log Z^r (which 
are centered) are finite and both asymptotic to ^ logiV as — > oo. 

Now we state the central limit theorem we shall apply (we follow page 87 
of Strook [19j). We assume that (X„) is a sequence of independent, square 
integrable valued random variables defined on the same probability space. 
Further we will assume that X„ has mean and strictly positive covariance 
cov {Xn). Finally we set: 

n n 
Sn = ^ Xm, Cn := COV {Sn) = ^ COV (X„) 
m=l m=l 

and 

S„ = (det(C„))i and5„ = -^. 

Theorem 4.2 (Multidimensional Central Limit Theorem, |19j p. 88). As- 
sume that 

A := lim 



exists and that 



1 " 

1™ ^ E ^ [\Xm\%X^\>e^„] = (4.8) 



n— »oo Yj, 

" m=l 



for every e > 0. Then the vector Sn converges in law to a Gaussian vector 
with mean and covariance matrix A. 
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Now we can prove Theorem 4.1 : 



N 



law 



ilogiV 



M + i7V2. 



Indeed, let us consider the Lyapounov exponents associated with (T„ 
and (Wn): 



L 



1 



N 



N 



and 



where sj^ 



L' 



1 



N 



N 



a 



3/2 

N n=l 



N 



E 



T->2 



and aj^ 



E 



3 



From the expressions of 



the cumulants, we have: 



,2 



2 



I ^ 1 



It is not hard to see, using the expression for the cumulants or the den- 
sity of the beta variables and the Wj, that the series E and 
X^^j^E [|W„p] both converge. Hence Lj^ and L'jy ^ as ^ cxd, 
and consequently (14. 8[) holds and the result follows from an application of 
the Multidimensional Central Limit Theorem of Theorem 14.21 

5. Iterated logarithm law 

In this Section, we give some iterated logarithm law for both the real and 
imaginary parts of the characteristic polynomial. Again, this can be done 
due to the decomposition given in Theorem 13. 2[ 

We first need some information about the rate of convergence in the cen- 
tral limit theorem. 



5.1. Further results about the rate of convergence. With the repre- 
sentation in Theorem 13.21 it is possible to obtain uniform and non-uniform 
estimates on the rate of convergence in the central limit theorem, using the 
Berry-Essen's inequalities (see [15] or [?]): 

Theorem 5.1. Let Xi, . . . be independent random variables such that 



E[Xj] = 0, andE[\Xj\ 



-1/2 



En 



^Xj<X 



'] < oo. Put (t| 
and 



E 



N 
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Then there exist two constants A and C not depending on n such that the 
following uniform and non uniform estimates hold: 

sup\Fn{x)-<^{x)\<ALn (5.1) 



and 



|F„ (x) - $ (x) I < 



CLr, 



(l + |x|) 



3 • 



(5.2) 



Now, applying the above theorem to the variables (Tj) and {Wj) we ob- 
tain: 

Proposition 5.2. The following estimates on the rate of convergence in the 
central limit theorem for the real and imaginary parts of the characteristic 
polynomial hold: 



log Zm 



logiV 



Dim log Z 



N 



< X 



< X 



loffiV 



$ (x) 



< 



< 



C 



i\ogNf^il + \x\f 
C 

(logiV)3/2(l + |x|)3 



(5.3) 
(5.4) 



where C is a constant. 

5.2. An iterated logarithm law. We first state a theorem of Petrov, [13^ 



Theorem 5.3 (Petrov). Let Xi, X2, ■ ■ ■ be independent random variables 
such that E [Xj] = 0, and a'j 



E 



^1 



< 00. Set Bn = Y.j=l CTj, Fn [x] 
Bn "^^^ Y^^=i -^j — ^ '^'^'^ ^{x) = J-00 c~*^^^dt. If the conditions 



1) Bn ^ 00 

Bn + l 



(2) 



Bn 



1; 



(3) sup, \Fn{x) - $(x)| = O ((log5„,) 
are satisfied for some 6 > 0, then 



-l-<5 



lim sup ■ 



1 a.s. 



(5.5) 



V2S„loglog5, 

Remark. If the conditions of the theorem are satisfied, then we also have: 

Sn 



lim inf ■ 



-1 a.s. 



^/2Bn log log B„ 

Before using Theorem 15.31 for the real and imaginary parts of logZ^v, 
we need to give the explicit meaning of the "almost sure convergence" for 
matrices with different sizes. 

hnagine that in Proposition 12.11 we choose for M the symmetry which 
transforms ei (the first vector of the basis) into Mi, a random vector of 
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^q'^^- Consider the set O = x x . . . endowed with the mea- 
sure ui X 1/2 X 1^3 ... , where is the uniform measure on the sphere 
(this can be a probabihty measure by defining the measure of a set as the 
hmiting measure of the finite-dimensional cylinders). Consider the applica- 
tion / which transforms uj £ O into an element of U{1) x U{2) x U{3) . . . 
with successive iterations of the Proposition 12.11 Then $7 = 3m{f) is nat- 
urally endowed with a probability measure = Jm(z^), and the marginal 
distribution of iJ-u on the k^^ coordinate is the Haar measures on U{k). 

Let (7 be a function of a unitary matrix U , no matter the size of U (e.g. 
g = det{Id—U)). The introduction of the set Q with measure nu allows us to 
define the "almost sure" convergence of {g{Uk),k > 0), where {Uk)k>o £ ^■ 
This is, for instance, the sense of the "a.s" in the following iterated logarithm 
law. 

Proposition 5.4. The following almost sure convergence (defined previ- 
ously) holds : 

hmsup— = = = 1, (5.6) 

V log A* log log log A* 

'3m\og Zn ^ 
limsup = = 1. (5-7) 

V log N log log log N 

Remark. The representations in law as sums of independent random vari- 
ables we have obtained could as well be used to obtain all sorts of refined 
large and moderate deviations estimates for the characteristic polynomial. 



6. Same results in the orthogonal case. 

The Mellin Fourier transform for Z := det{lN — M) found in [7] and [S] by 
Keating and Snaith, using the Selberg integrals, are, for the unitary group 
U{N) 

E nz|*e--*5^1 = fr T{k)T{k + t) 

and for the special orthogonal group S0{2N) 

E = 2^^* n —, ^- — i ^. (6.2 

L J iir(A:-i)r(t + A; + iV-l) ^ ^ 

Formula (j6.ip was directly proven in section [2j Here we show that such a 
probabilistic proof still holds for formula (|6.2p . 
Let 

^K^:={(ri,...,rjv)GM^: |ri|2 + ... + |rjvP = l} 

and fJ-o(N) be the Haar measure on 0{N). Then in analogy to Proposi- 
tion [2lTl we have: 
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Proposition 6.1. Let M G 0{N+1) he chosen such that its first column Mi 
is uniformly distributed on ,y^^^ . Let £ 0{N) be chosen independently 
of M according to the Haar measure ^o{N)- Then the matrix 

is distributed with the Haar measure /^o(iV+i)- 

The proof is essentially the same as the proof for Proposition 12. li If 
we choose for M a symmetry transforming ei in a uniformly chosen vector 
of we transform a random element of SO{n) into an element of 

0(n + 1) with determinant —1, and reciprocally. As a consequence, the 
following result, analogous of Proposition 12.21 can easily be shown. 

Corollary 6.2. Let SO £ S0{2n) be distributed with the Haar measure 
l^S0(2n)- Then 




with ei, . . . , e2n; /3i/2,i/2) • • • > /^i/2,(2n-i)/2 independent random variables such 
that P(efc = 1) = P(efc = —1) = 1/2, and the j3's being beta distributed with 
the indicated parameters. 

Remark. A direct calculation with the suitable change of variables shows 
that 

1 + Cfc^ / /?! fc-i = 2/3 fc-i fc-i , 

V 2 ' 2 2 ' 2 

from which formula (j6.2p can be easily recovered. 

Remark. The same reasoning can be applied to many other groups such 

as EI(n), the set of n x n matrices Hn on the field of quaternions with 
T 

Hn Hn = In- The symplectic group requires some additional work, to 
appear in a future paper. 
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